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Super resolution microscopy such as STORM and (F)PALM is now a well known method for biological 
studies at the nanometer scale. However, conventional imaging schemes based on sparse activation of 
photo-switchable fluorescent probes have inherently slow temporal resolution which is a serious limitation 
when investigating live-cell dynamics. Here, we present an algorithm for high-density super-resolution 
microscopy which combines a sparsity-promoting formulation with a Taylor series approximation of the 
PSF. Our algorithm is designed to provide unbiased localization on continuous space and high recall rates 
for high-density imaging, and to have orders-of-magnitude shorter run times compared to previous 
high- density algorithms. We validated our algorithm on both simulated and experimental data, and 
demonstrated live-cell imaging with temporal resolution of 2.5 seconds by recovering fast ER dynamics. 

Single-molecule localization microscopy methods, such as STORM 1 and (F)PALM 23 , utilize sparse activa- 
tion of photo-switchable fluorescent probes in both temporal and spatial domains. Each activated probe can 
be assimilated to an ideal point source so that the acquired images consist of isolated replicates of the point 
spread function of the microscope (PSF). This allows one to achieve sub-pixel accuracy on the order of tens of 
nanometers for the estimated location of each probe 1 " 6 . In general, reconstruction of sub -cellular structures relies 
on numerous localized probes, and the required acquisition time of these methods is therefore relatively long, i.e. 
on the order of minutes. This is a serious limitation when investigating live-cell dynamics. 

One possible approach for overcoming this limitation is high-density imaging 7 . By increasing the density of 
activated probes, shorter acquisition times for a single super-resolution image can be achieved. However this 
complicates the localization task due to overlapping PSFs. DAOSTORM 7 , for example, fits multiple overlapping 
PSFs in an iterative manner by analyzing pixel clusters in the residual image. The positions of the probes are 
determined by minimizing a least-squares criterion. CSSTORM 8 (Compressed sensing STORM) and 
deconSTORM 9 (deconvolution STORM) impose sparsity priors on the distribution of probes. In CSSTORM 
algorithm, the localization task is formulated as a convex optimization problem and solved by means of linear 
programming, while deconSTORM uses a modified Lucy- Richardson deconvolution algorithm by exploiting 
temporal correlation of activated probes. In general, these sparsity promoting methods provide increased recall 
rates compared to multi- emitter fitting at the expense of higher computational complexity. In a different 
approach, super resolution optical fluctuations imaging (SOFI) 10 and 3B analysis 11 utilize stochastic photon- 
emission processes such as photo-bleaching and blinking to reconstruct high-density data. For example, 3B 
analysis based on realistic models of photo-bleaching and blinking processes reconstructs the high-density data 
using a Bayesian approach. The Current implementation of 3B analysis has a relatively long reconstruction time 
but it can be made faster by using a computationally efficient Bayesian algorithm or parallel computing 12 . 

In addition, all of these sparsity-promoting methods are based on similar discrete formulations. They recon- 
struct a high resolution image on a pre-defined sub-pixel grid, e.g., with a pixel size of 20 nm. Such formulations, 
however, have three inherent limitations. First, discrete-domain formulations can account for only pre-defined 
locations, not all possible probe locations over a continuum. Therefore, these methods need to extract the 
localization information from their reconstructed high resolution images. A partial solution to this problem is 
to compute local centers of mass in the reconstructed image 8 , but this tends to result in a biased estimation of the 
probe locations. Second, using a finer sub-pixel grid increases the computational load, especially with the linear- 
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Algorithm schematic diagram 



1. Deconvolution with sparsity priors 

A deblurred image is generated by using sparsity-promoting priors (weighted 11 norm) 
on a sub-pixel grid. 




Deconvolution by using 
the weighted 11 penalty 



Take a spatial support of 
the deblurred image 



2. Deconvolution with fixed spatial support 

Deconvolution by minimizing the least-squares criterion on a fixed spatial support 



Deconvolution on the 
spatial support. 
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initial localization 
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3. Continuous refinement 

Initial localizations are refined by alternatively updating positions and brightnesses. 
Yellow arrows show directions of refinement. 
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Figure 1 | Schematic illustration for FALCON. For a detailed description, see Supplementary note. 



programming approach used in CSSTORM. The final issue is the 
underestimation of photon counts. Indeed, many sparse solvers are 
based on iterative thresholding, which often truncates signals toward 
zero, resulting in a spatial bias of the recovered locations of the probes. 

To address these limitations of current high-density localization 
algorithms, we introduce a FAst Localization algorithm based on a 
CONtinuous- space formulation (FALCON) for high density super- 
resolution microscopy data. In particular, to obtain a grid-free recon- 
struction, our approach combines a sparsity-promoting formulation 
with a Taylor approximation of the PSF. It consists of three main 
stages: deconvolution with a sparsity prior, deconvolution with a 
fixed spatial support, and continuous -domain refinement (Fig. 1 
and 2, and Supplementary Fig. 1, Supplementary note). The first step 
generates a high- resolution sparse image similarly to CSSTORM. 
The second deconvolution step corrects underestimated pixel values. 
The last stage uses a Taylor series approximation of the PSF for 



refining the position of each probe over a continuum. It also refines 
the photon count of each probe. For the Taylor expansions of the PSF 
with respect to the spatial offset from the grid location, we only use 
linear terms to make the optimization procedure faster. Importantly, 
this approximation is not limited to Gaussian functions. The final 
output of the algorithm are positions and photon counts of the 
localized probes. All these steps are implemented in a computation- 
ally efficient way by utilizing ADMM (Alternating Direction Method 
of Multipliers) 13 for fast reconstruction of the high resolution sparse 
images, and by using alternating refinements. 

Results 

We validated the performance of the proposed algorithm using both 
simulated and real STORM/PALM data. First, we confirmed that the 
linear approximation of the PSF provides a good approximation for a 
PSF that is shifted by a few nanometers from the sub -pixel grid. 
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Figure 2 | Taylor approximation of the shifted PSF for continuous refinement of FALCON, (a) First-order Taylor approximation of the shifted PSF: the 
shifted PSF is approximated by the PSF centered on the sub-pixel and the derivative of the PSF. (b) Approximation error of Gaussian PSF of 360 nm 
FWHM. The error of the approximation does not exceed 2% up to 30 nm displacement. 



Specifically, the relative approximation error is less than 2% for a 
displacement of up to 30 nm for a Gaussian function with a FWHM 
(Full Width at Half Maximum of 360 nm (Fig. 2b). Next, we inves- 
tigated the localization performance of our algorithm for a single 
probe. Starting from a uniform probability distribution over an area 



of a single sub-pixel of 33.3 nm (Fig. 3a-c), we localized 50, 000 
individual probes using two methods: center-of-mass on a decon- 
volved image 8 and our continuous refinement by Taylor series 
approximation. The former method resulted in estimations that were 
strongly biased towards the sub-pixel grid points, while our approach 
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Figure 3 | Numerical analysis of the continuous refinement on 50,000 simulated images. In each image, a single molecule is randomly placed with 
uniform probability distribution within a single sub-pixel area (a-c) or along a diagonal line (d-f). Histograms of ground-truth positions (a,d), 
histograms of the centroid fit that is applied to the high resolution images reconstructed by a sparsity-based deconvolution (U minimization) (b,e) and 
histograms of FALCON (c,f). 
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Figure 4 | Performances of FALCON in comparison with Least-square fitting, DAOSTORM and CSSTORM over a wide range of imaging densities. 

Simulation on the random distribution of molecules over a wide range of imaging densities with high-photon emission rates: recall rates (a), localization 
accuracy (b) and low-photon emission rates: recall rates (c), localization accuracy (d). The error bars represent standard deviations. 



produced a much more uniform distribution of the estimated particle 
locations. Similar results were also observed for probes that were 
uniformly distributed along a diagonal line (Fig. 3d-f). These simu- 
lations were performed with high photon-emission rates (5, 000 
photons on average, standard deviation of 2, 000). 

We also compared our algorithm with a least- square fitting 
method, DAOSTORM and CSSTORM over a wide range of imaging 
densities. We simulated two different photon-emission rates: high (5, 
000 photons on average, standard deviation of 2, 000) and low (350 
photons on average, standard deviation of 70). The FWHM of the 
PSF was fixed at 385 nm. For the high photon- emission rate simula- 
tions (Fig. 4a and b), FALCON and CSSTORM achieved two times 
higher recall rates than DAOSTORM did for densities larger than 
6 fim~ 2 . In terms of localization accuracy, the proposed algorithm 
outperforms all the other algorithms in high imaging densities (4- 
9.3 /mi~ 2 ), and our method is compatible with the least-square fit- 
ting showing only 5% difference in the lowest imaging density. For 
the low photon- emission rate simulations (Fig. 4c and d), the 
improvement is even more distinct in terms of recall rate, and 
FALCON demonstrates better localization accuracy over the whole 
range of imaging densities (0.05-6.5 /nm~ 2 ). 

The performance dependency on the size of PSF was also investi- 
gated. Fixing the imaging density to 6 ^m" 2 , we investigated the 
performance of FALCON, DAOSTORM and CSSTORM for various 
FWHMs in the range of 250-430 nm (Fig. 5). As the FWHM 
becomes smaller, the performance of all algorithms improves, 
leading to comparable results. Nevertheless, as the FWHM becomes 
larger, FALCON performs better than CSSTORM, leading to better 
accuracy and higher recall rates with a notable improvement in the 
low photon-emission cases (Fig. 5c and d). DAOSTORM always 
shows minimum recall rates and the differences with the others are 
very clear with the broader PSFs. 

In order to quantify the performance of our algorithm on specific 
geometric structures, we simulated circular structures whose radii 



varied between 150 nm and 400 nm (Fig. 6). In every frame, 4 acti- 
vated probes are randomly placed on the circle, resulting in overlap- 
ping PSFs even at the largest radius. We estimated the radius of the 
reconstructed circle by averaging the distances between all localized 
probes from 300 simulated images and the center of the circle. The 
results of the simulation demonstrate that our algorithm performs 
better for extracting this geometrically relevant measure, over the 
entire range of simulations. In contrast, the values by CSSTORM 
were underestimated by up to 18 nm as a result of underestimated 
photons; CSSTORM also shows spurious peaks as artifacts in the 
histogram (Fig. 6b) because of the discrete formulation as shown 
in (Fig. 3). These artifacts were not created by FALCON and 
DAOSTORM. In terms of the estimation bias, DAOSTORM was 
comparable to FALCON, but it had smaller recall rates and larger 
standard deviation of localized locations (Fig. 6c). 

Our ADMM implementation of the two deconvolution stages exhi- 
bits fast convergence, and the refinement stage is computationally 
efficient thanks to the closed-form expression of each update step. 
For example, on an Intel i7 3.4 GHz CPU, our stand-alone Matlab 
implementation of FALCON took 5 minutes to reconstruct 100 high- 
density (6 ^m" 2 ) STORM images of 100 ^m 2 (100 X 100 pixels); on a 
GPU (Nvidia GTX Titan) it took FALCON 40 seconds to reconstruct 
the same data. Since the run time of our method mainly depends on 
the raw camera size and over-sampling factor M for the reconstruction 
grid, the computation cost will be further reduced by using a coarser 
sub-pixel grid for the deconvolution steps. In FALCON, the role of the 
deconvolution steps is mainly to estimate the number of sources and 
initial locations, and these initial values can be refined. We demon- 
strated that FALCON using coarser grid of 50 nm provides compat- 
ible results at the small expense of recall rates in comparison with 
FALCON using the finer grid of 33.3 nm in spite of only using reduced 
number of variables to estimate as shown in Supplementary Fig. 5. 

We demonstrate the performance of FALCON on a real STORM 
data of fixed microtubules (MTs). FALCON reconstructed alpha 
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Figure 5 | Performances of FALCON in comparison with DAOSTORM and CSSTORM over a wide range of PSF widths from 250 nm to 430 nm. 

Simulated data with high-photon emission rates: (a) molecular recall rates, (b) localization accuracy and low-photon emission rates: (c) molecular recall 
rates, (d) localization accuracy. The error bars indicate standard deviations. 
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Figure 6 | Performance analysis for synthetic "ring" phantoms with various radii, (a) Rendered images of the algorithms overlaid by a cyan circle with 
150 nm radius as a ground truth, (b) Histogram of the distances from localized molecules by FALCON, DAOSTORM, and CSSTORM to the center of 
the ring in (a), (c) The differences between the true radius and the distances on average along various radii from 150 nm to 400 nm with 4 fixed molecules 
on the circle in every simulated image. The error bars represent standard deviations. 
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Figure 7 | FALCON performance on fixed microtubules data, (a) FALCON reconstructed a-tubulin subunits of microtubules (MTs) labeled with 
Alexa 647 in Cos-7 cells; 500 raw frames were used to reconstruct a super-resolution image. White box highlights the intersection of MTs reconstructed by 
FALCON, (b) Fourier ring correlation(FRC) analysis on the reconstructed images by DAOSTORM, CSSTORM and FALCON. For the FRC analysis, 
fiducial beads indicated by yellow circles are discarded. The error bars represent standard deviations. Cross-sectional profiles across region in white box 
(yellow line) were measured by DAOSTORM, CSSTORM and FALCON, and plotted in (c). Scale bars are 1 /urn in (a). 



tubulin subunits of microtubules labeled with Alexa 647 in Cos-7 
cells, by using 500 frames imaged at 25 Hz (Fig. 7a). For the quant- 
itative analysis, as increasing a number of raw frames to reconstruct a 
image, we measured resolutions of the reconstructed images by 
means of Fourier ring correlation (FRC) criteria 14 . For the FRC ana- 
lysis, each half of localized locations were binned into an image of 
5nm pitch, and cross correlation between the two images was calcu- 
lated, then converted as a FRC resolution. FALCON shows a reso- 
lution improvement of up to 20% over the other high-density 
algorithms (Fig. 7b). In other words, FALCON can achieve the same 
spatial resolution using 100 times fewer frames. The superiority of 
FALCON is especially visible at the intersection of MTs. FALCON 
offers an unbiased estimation of distances between MTs, showing 
fewer spurious and noisy peaks than other high density algorithms 
(Fig. 7c). This is in accordance with our results concerning the radii 
of simulated circles. FALCON also provides a reliable distribution of 
photon counts which are well matched to results of another low- 
density MT data, acquired under identical preparation protocols 15 , 
whereas DAOSTORM (green line) loses many molecules with 
photon counts < 5,000 in comparison and CSSTORM yields many 
false-positive molecules with photon counts <500 (Supplementary 
Fig. 12). Furthermore, the proposed method showed improvements 
even in low- density data; especially around regions of complex MT 
structures or inserted fiducial beads, in comparison with Least- 
square fitting method (Supplementary Fig. 13). 

We also applied FALCON to reconstruct live, PALM data. In par- 
ticular, the dynamics of the endoplasmic reticulum (ER) were captured 
in U20S cells labeled with tdEos fused to reticulon-4; reconstructed 
movies were made, with each frame being a single PALM image 
acquired over 2.5 seconds (160 frames of raw data). Over the course 
of 20 seconds, we observed several changing features of the ER net- 
work. In particular, expanding and shrinking tubules were visualized 



along with the dynamic motion of tubule junction points (Fig. 8f-j, 
white arrow). We also captured distinct features of this organelle 
reshaping itself: the disappearance of tubule junctions (Fig. 8f-j, dotted 
circle) and emergence of tubules (Fig. 8f-j, blue arrow) to form new 
junction points are just a few examples of this. Importantly, we mea- 
sured an average tubule thickness of approximately 60 nm (Fig. 8g), 
which is in line with previously reported diameter measurements 1617 . 
Furthermore, we compared super-resolution images reconstructed by 
DAOSTORM, CSSTORM and FALCON by increasing accumulation 
time scales (1-3 s) in Supplementary Fig. 14. FALCON shows better 
reconstruction of complex ER structures than the others, and detected 
more molecules, up to 50% in comparison with DAOSTORM and 
CSSTORM, which can improve temporal resolution. 

Discussion 

We presented a fast localization algorithm for high-density data which 
incorporates a sparsity promoting formulation together with a Taylor 
series approximation of the PSF for continuous localization. The 
robustness of the proposed method has been extensively studied by 
simulated and experimental data in comparison with previous high- 
density algorithms based on multi- emitter fitting (DAOSTORM) or 
deconvolution using sparsity-promoting priors (CSSTORM and 
deconSTORM). Localization performances have been quantitatively 
analyzed not only in terms of recall rates and localization accuracy, 
but also in terms of localization bias. Importantly, even though local- 
ization bias can be easily observed in case of high-density imaging and 
is important for quantitative biological analyses, it has not been prop- 
erly investigated before 7 " 9 . 

We have shown that a localization bias can be induced by several 
factors. Low recall rates is one factor, which can occur by counting 
multiple closely-spaced probes as one, and then localizing the cen- 
troid position. In general, a deconvolution based approaches utilizing 
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Figure 8 | FALCON performance on live ER data. Live imaging of the endoplasmic reticulum protein, reticulon-4 fused to tdEos imaged over 20 seconds 
in a U20S cell. Conventional (a-e) and super-resolution (f-j) snapshots are shown with a 2.5 seconds temporal resolution. The average size of tubules 
measured from the reconstruction is approximately 60 nm (at FWHM). Yellow markers in (g) highlight a representative tubule width. The dynamic 
motions of these structures are highlighted in (g-j) indicated by white and blue arrows and dashed, yellow circles. All scale bars are 1 /urn. 



a sparsity promoting priors can achieve higher recall rates than 
multi- emitter fitting based methods by effectively resolving closely 
placed probes. However, this approach can introduce an additional 
bias for location estimation mainly originating from its use of the 
discrete sparsity-promoting formulation. Our experiments have illu- 
strated grid-pattern artifacts and spatial distortions appearing in the 
final image reconstructed when using the discrete sparsity-promot- 
ing formulation. In particular, spatial distortions can be problematic 
when imaging intracellular structures. For instance, quantifying the 
sizes of mitochondria, microtubules or the tubules of the ER would be 
done with less confidence, since structures would be falsely con- 
stricted. On the other hand, the proposed method based on the 
continuous formulation minimized these bias errors and has shown 
almost bias-free localization results. 

In live imaging, the variance of the localization accuracy is com- 
posed of both the algorithmic localization variance and the variance 
of the fluctuations. Indeed, the motion artifacts due to dynamic 
structures that move on time-scales shorter than our exposure time 
can be problematic, often resulting in artificially larger or smaller 
structures- as previously reported 18 . In the ER experiment, we sought 
to probe the dynamic motion of ER tubules. Previous studies 19 ' 20 , 
with similar exposure times demonstrated that motion of this organ- 
elle is dynamic on longer timescales than the one used here. 

Our method is currently using a two dimensional PSF model 
which is appropriate to a relatively flat sample region having a thick- 
ness below 1 /urn. However, it must be noted that our continuous 
sparsity-promoting formulation is not limited to a two-dimensional 
Gaussian PSF model. It can be easily extended to non-Gaussian 
models including three dimensional PSF models. For example, 
multi-plane imaging systems with three-dimensional PSF engineer- 
ing techniques can be accommodated with our formulation. 

The proposed formulation with a Gaussian noise assumption can 
be further improved by considering more realistic noise statistics. 



Since many state-of-the-art cameras for localization microscopy, 
such as an EMCCD (Electron multiplying charge -coupled device) 
or sCMOS (scientific complementary metal-oxide-semiconductor) 
camera, have very small read-out noise 21 , Poissonian shot noise 
becomes a dominant noise source to resolve. Basically, exploiting 
Poisson statistics 22 would be more advantageous especially for live- 
imaging data, which is likely to be achieved at low SNR due to low 
emission photons and high background. However, even with our 
current implementation based on a Gaussian noise model, at the 
SNR levels from our experimental conditions, the proposed algo- 
rithm is robust under Poisson noise statistics. 

Concerning a modeling of readout noise for our simulations, we 
assumed that a CCD -type camera was used, where readout noise is 
known to be well approximated by the normal distribution. The well- 
established high-density algorithms such as DAOSTORM and 
CSSTORM also assumed Gaussian readout noise of small variance 
for their simulations 7 ' 8 . Although sCMOS cameras have readout 
noise with pixel-to-pixel variations resulting in a non-Gaussian dis- 
tribution such as the log- normal distribution 23 , our simulation 
results indicate that the proposed algorithm is not sensitive to the 
choice of the noise model (Supplementary Fig. 9). 

In general, it is difficult to make a fair comparison between algo- 
rithms in terms of run time due to the difference in implementation 
environment. Since the speed of the algorithm is determined by 
various factors including optimization, programming language, 
operating system and usage of parallel implementation, the imple- 
mentation environment must also be thoroughly evaluated to 
maximize calculation efficiency. For example, a computationally 
optimized CSSTORM implementation 24 is faster by two orders of 
magnitude than the initial MATLAB -based implementation 8 which 
takes approximately one day to reconstruct a 100 fim 2 region of 100 
raw frames of STORM data with high-density activation(6 /nm~ 2 ) on 
an Intel i7, 3.4 GHz CPU. 
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Although the speed of our approach is also dependent on the 
implementation environment, a generally fast algorithm speed can 
be achieved since the algorithm starts from a coarser grid which 
inherently reduces the problem complexity. Specifically, the pro- 
posed method would be very beneficial for three dimensional local- 
ization problems with high-density data, where calculation time 
reduction is critical since discrete formulation requires numerous 
finer voxels to be estimated. 

In conclusion, Taylor approximation of the PSF can be jointly 
utilized with sparsity-promoting formulation for high- density 
imaging to yield continuous localization with reduced complexity 
(up to several orders of magnitude). Furthermore, we have also 
investigated how a spatial bias is introduced in the localization 
estimation of high-density imaging. We also have showed that the 
proposed continuous formulation with specially designed sparsity 
priors substantially reduces the bias and provides better localization 
accuracy and higher recall rates than other currently available meth- 
ods over a wide variety of experimental conditions. Moreover, our 
method is efficient in terms of implementation, which reduces com- 
putation time by orders of magnitude in comparison with the pre- 
vious high density algorithms. These results were well confirmed 
through experimental data obtained using fixed microtubules sam- 
ple. We also successfully reconstructed fast dynamic motions of ER 
and measured its cross -sectional profiles reliably. 

Methods 

Simulation & Evaluation. 

Photon- emission statistics: We assume that the total number of photons emitted from 
a probe is a random variable that follows a log-normal distribution 8 . We considered 
two photon-emission rates: high and low. The former setting was primarily meant to 
simulate STORM data while the latter was closer to PALM data. The log-normal 
parameters for the high photon-emission case are (mean, standard deviation) = 
(5,000, 2,000), and for the low photon-emission case, they are (350, 70). Constant 
background fluorescence photons of 100 or 10 are added to every camera pixel for 
high and low emission cases, respectively. 

PSF model: When reconstructing real STORM data, we used low-density STORM 
images for fitting the measured PSF. Our PSF model is a sum of two Gaussian 
functions 8 . Low- density data is used as a reference of the PSF model, in order to 
optimize their widths and weights. For our simulated data, we chose to have two 
Gaussian functions that have a standard deviation of er and 2a, and that are linearly 
combined with a ratio of 4 : 1. 

Noise statistics: For simulating noise, we assumed that a high-performance CCD 
camera was used for data acquisition and considered two type of noise contribution, 
shot noise and a small readout noise 21 ' 23 . The shot noise follows Poisson statistics. 
Specifically, the mean of the detected photons is determined by the expected number 
of incoming photons and by the background fluorescence level at each pixel. We also 
added CCD -type readout noise which follows a normal distribution with zero mean 
and variance of 2. In Supplementary Fig. 9, we also considered a CMOS-type camera 
whose readout noise has intrinsic pixel to pixel variation, resulting in non- Gaussian 
distribution. Specifically, CMOS-type readout noise follows a log- normal distribution 
with unit mean and variance of 2. 

Probe distribution and performance analysis: For the experiment of Fig. 4 and 5, we 
randomly distributed N probes in the central region of the image; the size of the region 
is 100 X 100 pixels. The size of the image is slightly larger, 108 X 108 pixels. The 
camera pixel size is 100 nm. In Fig. 4, the imaging densities vary between 0.05 /mi" 2 
and 9.3 /mi" 2 for the high photon-emission rates (Fig. 4a and b) and between 
0.05 /mi" 2 and 6.5 /mi" 2 for the low photon-emission rates (Fig. 4c and d). The 
FWHM of the PSF was set to 385 nm. We used 50 realizations for every density. Fig. 5 
depicts additional analysis with FWHM PSF between 250 nm and 430 nm. The 
imaging density is 6 /mi" 2 and we used 30 realizations for each FWHM. 

The performance analysis was carried out by matching each localized probe with 
the closest ground-truth probe. Localization errors larger than 300 nm were excluded 
from the analysis. The localization accuracy is expressed in terms of standard devi- 
ation or the FWHM of the histograms of localization errors. The recall rate is defined 
by the ratio between the number of matched probes and N. 

Simulation of geometric structures: We generated STORM images by distributing 
particles randomly on a circle of radius r. The radius varies between 150 nm and 
400 nm. The probes are activated at different time instances, covering 300 frames. 
The number of probes at each frame is N = 4 (Fig. 6). Supplementary Fig. 4 depicts 
estimation results for various number of molecules N between 1-6 with fixed 
r = 200 nm. For every localized probe, we calculated the distance from the probe to 
the center of the circle, and estimated the radius by averaging all distances. Then, we 



plot the differences between the true radii and the estimated values. Here, high 
photon-emission rates are used. 

Image and signal processing. 

We used Gaussian rendering to generate super- resolution images. For a rendered 
image, every localized molecule is convolved with a Gaussian kernel, and then scaled 
by its estimated photon count. In Fig. 8, normalized Gaussian kernels were used 
without photon scaling. For the experimental STORM data of fixed microtubule 
(Fig. 7), stage drifts are corrected by using fiducial beads. To correct the drift of the 
sample in Supplementary Fig. 13, the mean vertical (resp. horizontal) position of a 
straight horizontal(resp. vertical) segment of microtubules was measured as a 
function of frame number, and the polynomial fit of this function was subtracted from 
the vertical (resp. horizontal) coordinates of all the peaks. 

In order to determine the achievable temporal resolution of our live data, we used 
the Nyquist criteria (p = (2 /R) ) as commonly applied to super-resolution imaging 25 . 
We counted the number of molecules on our structure of interest and obtained an 
estimate of its effective area by thresholding the super-resolution image; from this we 
had an estimate of our localization density. With a 2.5 s reconstructed super-reso- 
lution image, our localization density corresponds to R = 40 nm. 

Sample preparation & Data acquisition. 

Cell culture: COS-7 and U20S cells were cultured in DMEM supplemented with 10% 
FBS (Sigma Aldrich) in a cell culture incubator (37°C and 5% C02) and plated at low 
confluency on cleaned 25 mm size 1 cover-glass (Menzell). 

Sample preparation for microtubule imaging infixed cells: Prior to fixation, all solu- 
tions were pre-warmed at 37°C: 24 hours after plating, Cos-7 cells were pre-extracted 
for 10 s in 0.5% Triton X-100 (Triton) in BRB80 (80 mM PIPES, 1 mM MgCl 2 , 

1 mM EGTA, adjusted to pH 6.8 with KOH) supplemented with 4 mM EGTA, 
washed in PBS, fixed for 10 min in - 20°C-Methanol (Sigma), and washed again in 
PBS. The samples were then blocked 30 minutes in 5% BSA, before being incubated 
for 1.5 hour at room temperature with 1 : 1000 mouse a-tubulin antibodies (Sigma, 
T5168) in PBS - 1% BSA - 0.2% Triton (PBST), followed by 3 washes with PBS-0.2% 
Triton, and then incubated for 45 min in PBST with 1 : 1000 goat anti-mouse Alexa- 
647 F(ab)2 secondary antibody fragments (Life Technologies, A-21237). 

Sample preparation for live-cell imaging of the endoplasmic reticulum (ER): After 
letting cells grow to approximately 60% confluency, U20S cells were transfected with 

2 fig of the reticulon 4 (Rtn4)-tdEos plasmid using FuGENE 6 transfection reagent 
(Roche). Cells were washed 24 hours post transfection and imaged in Leibovitz. 

Live and fixed super-resolution imaging: Imaging of the live and fixed samples was 
performed on a modified Olympus 1X71 inverted microscope. In the case of Alexa- 
647 labeled a-tubulin, a laser at 641 nm (Coherent, CUBE 640- 100C) was reflected by 
a multi-band dichroic (89100 bs, Chroma) onto the back aperture of a 100 X 1.3 NA 
oil objective (Olympus, UplanFL) mounted on a piezo objective scanner (P-725 
PIFOC, Physicinstrument). The collected fluorescence from the sample was filtered 
by a band-pass emission filter (ET700/75, Chroma) and imaged onto an EMCCD 
camera (IxonEM+ , Andor) with a 100 nm pixel size and using the conventional CCD 
amplifier at a frame rate of 25 fps. Laser intensity on the sample was — 1-2 kW. cm" 2 . 
A 405 nm laser was added to maintain a high peak density. Blinking was induced 
according to 15 . Briefly, 30 pL of a mixture of 20% Vectashield, 80% Glycerol was 
pipeted on top of the sample, and a cover- glass was added on top of the buffer to 
spread it evenly on the sample. For live-cell imaging, transfected cells were identified 
using a 488 nm laser (Sapphire 488, Coherent) at an intensity < 1 kW cm" 2 to 
prevent activation of tdEos. Subsequently, cells were imaged using a 561 nm laser 
(Sapphire 561, Coherent) with an intensity of approximately 3 kW cm" 2 . A 405 nm 
laser was used to convert tdEos molecules from their green to red form and was left on 
throughout acquisition. Fluorescence was collected after passing through a ET605/70 
(Chroma) emission filter. Imaging was performed at a frame rate of 64 fps. 
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